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Abstract 

We analyze the observed properties of dwarf galaxies, which are dark matter dominated astro- 
physical objects, by assuming that dark matter is in the form of a strongly - coupled, dilute Bose 
- Einstein condensate. The basic astrophysical properties of the condensate (density profile, rota- 
tional velocity, and mass profile, respectively), are derived from a variational principle. To test the 
validity of the model we compare first the tangential velocity equation of the model with a sample 
of eight rotation curves of dwarf galaxies. We find a good agreement between the theoretically pre- 
dicted rotation curves (without any baryonic component) and the observational data. The mean 
value of the logarithmic inner slope of the mass density profile of dwarf galaxies is also obtained, 
and it is shown that the observed value of this parameter is in agreement with the theoretical 
results. The predictions of the Bose - Einstein condensate model are also systematically compared 
with the predictions of the standard Cold Dark Matter model. The non-singular density profiles of 
the Bose-Einstein condensed dark matter generally show the presence of an extended core, whose 
presence is due to the strong interaction between dark matter particles. 
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I. INTRODUCTION 

Despite the important achievements of the A Cold Dark Matter (ACDM) model at 
galactic scales ~ 10 kpc, the standard cosmological model meets with severe difficulties in 
explaining the observed distribution of the invisible matter around the luminous one. In 
fact, A^-body simulations, performed in this scenario, predict that bound halos surround- 
ing galaxies must have very characteristic density profiles, called the Navarro- Frenk- White 



(NFW) profiles, which feature a well pronounced central cusp 



a, 
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where r s is a scale radius and p s is a characteristic density. On the observational side, high- 
resolution rotation curves show, instead, that the actual distribution of dark matter is much 
shallower than the above, and it presents a nearly constant density core |3f] . 

The discrepancies between observations and simulations done in the framework of the 
ACDM model are particularly important in the case of the recent observations of the dwarf 



galaxies done by "The Hi Nearby Galaxy Survey" (THINGS) J4]. Dwarf galaxies are dark 
matter dominated cosmic structures, with a very small contribution of baryons to the total 
matter content. Similarly to the Low Surface Brightness (LSB) galaxies, they are ideal 
objects for the study of the dark matter properties. The observed rotation curves of the 
dwarf galaxies rise too slowly, as compared to the rotation curves derived from the cusp- 
like dark matter distribution in CDM halos. Instead, they are better described by core-like 
models, dominated by a central constant-density core |5j. Moreover, the mean value of the 
ogarithmic inner slopes of the mass density profiles of the dwarf galaxies is a = —0.29 ±0.07 
5], while for a large sample of LSB galaxies a = —0.2 ±0.2 [6[, values which are significantly 
different from the steep slope of a = —1.0, inferred from dark-matter-only simulations. A 
second serious problem is related to the mass of the dwarf galaxies obtained from simulations. 
The dwarf galaxies formed in current hydrodynamical simulations are more than an order 
of magnitude more luminous than expected for haloes of masses of the order of 1O 1O M Q. 
Therefore explaining the observed properties of dwarf galaxies has become a very challenging 
task, and "missing astrophysical effects in the simulations are the most likely cause of the 
discrepancy, and the most promising target in search of its resolution" Q. 

Recently, it was suggested that cold dark matter particles interacting through a Yukawa 
potential could produce the cores in dwarf galaxies [8j. Simulations that include the effect 



of baryonic feedback processes, such as gas cooling, star formation, cosmic UV background 
heating and physically motivated gas outflows driven by supernovae have succeeded in ob- 
taining a value of a = —0.4 ±0.1 [9]. 

The possibility that dark matter is in the form of a Bose-Einstein Condensate (BEC), 
namely, an assembly of light individual bosons that acquire a repulsive interaction by oc- 
cupying the same ground energy state, has been extensively investigated in the physical 



literature 
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J-llJ]. The experimental realization of the Bose-Einstein condensation in ter- 
restrial laboratories for a large class of particles (atoms) and physical systems has given a 
sound theoretical and experimental support to this hypothesis. 

The Bose-Einstein condensation process was first observed experimentally in 1995 in 
dilute alkali gases, such as vapors of rubidium and sodium, confined in a magnetic trap, and 
cooled to very low temperatures. A sharp peak in the velocity distribution was observed 
below a critical temperature, indicating that condensation has occurred, with the alkali 
atoms condensed in the same ground state and showing a narrow peak in the momentum 
space and in the coordinate space [l5|. 

Different properties of the Bose-Einstein condensed dark matter were investigated recently 



in 16 - 



271 ] . It was also shown that the creation of quantum vortices are favored in strongly- 



coupled condensates, while this is not the case for axion condensates, in which the particles 
are effectively non- interacting j^ . 

The possibility that the core-cusp problem of standard CDM models can be solved by 
assuming that dark matter is composed of ultralight scalar particles, with masses of the 
order of 10 -22 eV, initially in a cold Bose-Einstein condensate ( "fuzzy dark matter"), was 



considered in 



28] . The wave properties of the dark matter stabilize gravitational collapse, 



providing halo cores. In this model stability below the Jeans scale is guaranteed by the 
uncertainty principle, since an increase in momentum opposes any attempt to confine the 
particle further. The one- dimensional simulations performed in 28j suggest that the density 
profile on the Jeans scale is not universal, but also evolves continuously on the dynamical 
time scale (or faster), due to quantum interference effects. 

An analysis of the rotation curves in the Gross-Pitaevskii equation based BEC model 
of dark matter was performed, for several High Surface Brightness (HSB), LSB, and dwarf 
ga , X1 es, respective,, ta Q. There is a good general agreement between the theoretical 
predictions and observations. However, there are several reasons to extend the analysis 
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initiated in |l4j. First, there is a significant increase in the observational accuracy of rotation 
curve data for dwarf galaxies, which are dark matter dominated objects, and more high 
precision data are avai^e present* QQ. Second the, is a major discrepancy between 
observations and the theoretical predictions of the CDM models for the mean value of the 
logarithmic inner density slopes of the dwarf galaxies. The mean values of the logarithmic 
inner density and velocity slopes are observed quantities that could prove more useful in 
discriminating between different dark matter models and observations than the rotation 
curves themselves, which can be fitted by many models. Therefore it is important to analyze 
these quantities in the framework of the BEC models. Thirdly, the explicit presentation and 
detailed discussion of the BEC density profiles for a sample of dwarf galaxies has not been 
done yet. And finally, a systematic comparison between the prediction of the Navarro - 
Frenk - White dark matter density profiles, given by Eq. ([1]), and the BEC profiles, is still 
missing. 

It is the main goal of the present paper to address the above mentioned issues, by ana- 
lyzing the observed properties of the dark matter halos of eight dwarf galaxies. The dark 
matter is assumed to be in the form of a strongly - coupled, non-relativistic diluted Bose- 
Einstein condensate, as proposed in |14J. As a first step in our study we derive the basic 
astrophysical parameters of the condensate (density, velocity and mass profiles, respectively) 
from a variational principle. Then the theoretical predictions of the model are compared 
with the observational data for eight dwarf galaxies. The theoretical predictions of the ro- 
tation curves are compared with the observed rotation curves of the dwarf galaxies, as well 
as with the results derived from the NFW profile. There is a general good agreement be- 
tween the theoretical predictions of the BEC model and the behavior of the rotation curves. 
The observed mean value of the logarithmic inner slope of the mass density profile can also 
be theoretically reproduced. The dark matter density profiles for these galaxies are also 
explicitly presented, and compared with the standard NFW profiles. 

The present paper is organized as follows. In Section [III the basic equations of the BEC 
model are written down. The predictions of the BEC and NFW models are compared with 
the observations in Section II III We discuss and conclude our results in Section IIVI 
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II. BOSE-EINSTEIN CONDENSED DARK MATTER 



At very low temperatures, all particles in a dilute Bose gas condense to the same quantum 
ground state, forming a Bose-Einstein Condensate (BEC). Particles become correlated with 
each other when their wavelengths overlap, that is, the thermal wavelength A^ is greater 
than the mean inter-particles distance I. This happens at a temperature T < 2irh 2 /mkBn 2 ^ 3 , 



where m is the mass o: 
Boltzmann's constant 



the particle in the condensate, n is the number density, and k^ is 
291 ] . A coherent state develops when the particle density is enough 
high, or the temperature is sufficiently low. We assume that the dark matter halos are 
composed of a strongly - coupled dilute Bose-Einstein condensate at absolute zero. Hence 
almost all the dark matter particles are in the condensate. In a dilute and cold gas, only 
binary collisions at low energy are relevant, and these collisions are characterized by a single 
parameter, the s-wave scattering length a, independently of the details of the two-body 
potential. Therefore, one can replace the interaction potential with an effective interaction 
Vj (f — f) = X5 (f — r), where the coupling constant A is related to the scattering length a 
through A = Aiih 2 a/m [29]. The ground state properties of the dark matter are described 
by the mean-field Gross-Pitaevskii (GP) equation. The GP equation for the dark matter 
halos can be derived from the GP energy functional, 



E 



r |^|V^(r)| 2 + ^|^(r)| 4 
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where ip (r) is the wave function of the condensate, and Uq = <AnTi 2 a/m [29]. The first term 
in the energy functional is the quantum pressure, the second is the interaction energy, and 
the third is the gravitational potential energy. The mass density of the condensate is defined 
as p(f) = m \ip (r)| 2 , and the normalization condition is N = J \ip (f)| 2 df, where N is the 
total number of dark matter particles. The variational procedure 5E [ijj] —fi5 J \tp (r) | 2 df = 
gives the GP equation as 

h 2 

~ TT-^ ^ + V (r)4> if) + U \i){f)\ 2 i){f) = /xV> (f) , (3) 

where /i is the chemical potential, and the gravitational potential V satisfies the Poisson 
equation V 2 V = AnGp. When the number of particles in the condensate becomes large 
enough, the quantum pressure term makes a significant contribution only near the boundary, 



and is much smaller than the interaction energy term. Thus the quantum pressure term 
can be neglected (the Thomas- Fermi approximation). When iV — > oo, the Thomas- Fermi 
approximation becomes exact [29!]. Therefore we obtain 
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The Poisson equation becomes 



P if) = — [fi - mV (r)] . 



V 2 p + k 2 p = 0, 



(4) 



(5) 



where k = y AnGm 2 /U$. 

Hence the density distribution p of the static gravitationally bounded single component 
dark matter Bose-Einstein condensate is given by [14| 

sin kr 



P (r) = pc 



kr 



(6) 



where p c is the central density of the condensate, p c = p(0). The density profile is non- 
singular at the center r = 0. The mass profile m(r) = 4tt Jq p(r)r 2 dr of the galactic dark 
matter halo is 

zLtt nr 

(7) 



Airpr 

m (r) = (1 — kr cot kr) 



k 2 

with a boundary radius R. At the boundary of the dark matter distribution p(R) = 0, 
giving the condition kR = n, which fixes the radius of the condensate dark matter halo as 



R — 7T1 



' h 2 a 
Gm 3 



(8) 



The tangential velocity V of a test particle moving in the condensed dark halo can be 
represented as 



V^(r) = ^1 = ^ (I -kr cot kr) 



or, alternatively, 



V (r) 



AGp c R 2 
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k 2 

sin {nr / R) 
nr/R 



cos 



nr 



The total mass of the condensate dark matter halo M can be obtained as 

3 / 2 a 



M = 47T' 



( h 2 a\ 

\G^) P <=n Rp c 



(9) 
(10) 

(11) 



giving for the mean value < p > of the condensate density < p >= 3p c /n 2 . Hence the 
tangential velocity can also be represented as 



V (r) 



GM 
~R~ 



sin (nr / R) 
nr/R 



cos 



nr\ 
~RJ 



(12) 
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For r > R the rotation curves follow the standard Keplerian law. The mass of the particle 
in the condensate is given by 

(2fc2 \ 1/3 
I w 6.73 x 10^ 2 [a (fm)] 1/3 [R (kpc)]^ 2/3 eV. (13) 

For a « 1 fm and i? ~ 10 kpc, the typical mass of the condensate particle is of the order 
of m ~ 14 meV. For a ~ 10 6 fm, corresponding to the values of a observed in terrestrial 
laboratory experiments, m ~ 1.44 eV. These values are consistent with the limit m < 1.87 



eV obtained for the mass of the condensate particle from cosmological considerations 



30]. 



III. PROPERTIES OF DWARF GALAXIES IN THE BEC MODEL 



In order to analyze the properties of the dwarf galaxies in the framework of the BEC and 
CDM models we have selected a sample of eight dwarf galaxies for which high resolution 
rotation curves are available. The galaxies we consider are IC 2574, NGC 2366, Holmberg I 
(Ho I), Holmberg II (Ho II), M81 dwB, DDO 53 Q, DDO 39 [31 1 and DDO 154 Q, respec- 
tively. As a first step in the comparison of the observed properties of dwarf galaxies with 
the BEC model and the CDM model, respectively, we analyze the behavior of the rotation 
curves, by comparing the predictions of Eq. (jH]) and of Eq. ([TJ) with the observed rotation 
curves of the eight galaxies. This comparison allows us to determine either theoretically 
(BEC model), or by fitting (NFW model) the free parameters in the equations. Then, as a 
second step in our study we obtain the dark matter density profiles of the galaxies, and we 
estimate the mean values of the logarithmic density and velocity slopes in the BEC model. 



A. Rotational velocity profiles of the dwarf galaxies 

We begin our study of the BEC dark matter halo model by comparing its theoretical 
predictions with the observational data of the rotation curves of eight dwarf galaxies, as well 
as with the velocity curves obtained by using the standard NFW profile, given by Eq. ([Q). 
For the NFW density profile the rotational velocity is given by 
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Galaxy 


IC 2574 


NGC 2366 


Ho I 


Ho II 


DDO 53 


M81 dwB 


DDO 39 


DDO 154 


R (kpc) 


12.6 


6.5 


2.2 


7.0 


2.54 


1.0 


9.5 


7.5 


p c (lO- 24 g/cm 3 ) 


0.44 


0.89 


3.7 


0.35 


1.8 


18.0 


1.04 


0.60 



TABLE I: Radii and central densities of the dwarf galaxies used in the predicted BEC rotation 
curves. 



Galaxy 


IC 2574 


NGC 2366 


Ho I 


Ho II 


DDO 53 


M81 dwB 


DDO 39 


DDO 154 


p s ( g/cm 3 ) 


25272.5 


48553.7 


45877.9 


32341.2 


38835.4 


73763.8 


57633.0 


37276.1 


r s (km) 


49407.1 


29113.2 


27598.2 


18885.3 


3030.7 


47856.5 


35439.0 


22927.1 



TABLE II: The numerical values of the NFW parameters p s and r s obtained by fitting the NFW 
rotational velocity profile. 

In order to obtain the rotation curves in the BEC model we use the following procedure. 
The BEC theoretical rotation curves are not obtained by fitting the observational data. In- 
stead, a reasonable value for the radius R, interpreted as the point where the rotation curve 
starts to decline, is obtained from observations, and then a central density value is assumed. 
Therefore the theoretical rotation curves are predictions of the model, and not fits. On the 
other hand, we fit the NFW profiles with the observational data. 

The rotation curves predicted by the BEC model and the NFW profile, respectively, are 
presented, together with the observational data, in Fig. [TJ The numerical values of the two 
free parameters of the BEC model (central density and radius) are shown in Table [H 

The parameters obtained by fitting the NFW model with the observational results are 
presented in Table ILT1 

As one can see from Fig. (TJ there is an overall good agreement between the observed data 
and the theoretical predictions of the BEC model. The good agreement with the observations 
strongly suggest that our model may be relevant for obtaining a correct description of the 
dark matter, and its properties. In all predictions and fits the effect of the baryonic matter 
was completely neglected. Even that baryonic matter may represent only 10% of the entire 
mass of the dwarf galaxy, by including its gravitational effects, one can obtain a significantly 
better description of the rotation curves, especially in the small r regions. One should also 
mention that the model gives a relatively good description of the rotation curve of the galaxy 
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DDO 154, which is generally difficult to be achieved. However, the considered cases shows a 
variety of behaviors. Below we discuss the results obtained for the rotation curves for each 
galaxy. 

• IC 2574: For this galaxy the BEC model gives an almost perfect prediction for the 
rotation curve. The NFW velocity also gives a good fit. However, is is not as good as the 
BEC model prediction. The galaxy may be affected by non-circular motions |5|, but the 
rotation velocity, with a maximum of around 80 km/s is the highest in the sample. The 



galaxy is also the most massive in the considered sample, with M = 1.462 x 10 M ( 



© 



fl. 



With the use of the BEC model mass - radius relation, M = (4/7r) p c R 3 , we obtain for the 
mass of the galaxy the value M^^g 74 = 1.64 x 1O 1O M , which is very close to the observed 
value. 

• NGC 2366: The rotation curve of this galaxy is rather complicated, indicating that 
disturbances caused by non-circular motions could be present in the outer regions with 
r > 5 kpc. Non-rotational motions may also be present for r < 5 kpc, distorting the inner 
region of the rotation curve. It is difficult to give an estimate of the radius of the dark 
matter distribution, because of the increasing trend in the velocity field for r between 6 and 
8 kpc. The adopted values for the radius and central density give an estimate of the mass 
of the galaxy as M^g 2365 = 4.54 x 1O 9 M , while the observed mass is M = 4.29 x 1O 9 M 
jlj]. The BEC velocity profile gives a good description of the rotation curve for r 6 (4.5, 6.5) 
kpc. The NFW fit gives a reasonable good description of the rotation curve. 

• Holmberg I: For r > 1.5 kpc, the rotation curve of the galaxy Holmberg I shows a severe 
distortion, probably caused by the presence of strong external gravitational interactions, 
generating non-circular motions. Therefore the physical parameters and behavior of the 
galaxy without the presence of external perturbation are difficult to be extracted from the 
observed data. With the adopted values for the radius and central density, R = 2.2 kpc 
and p c = 3.7 x 10 -24 g/cm 3 , the mass of the galaxy is obtained as M^pf = 0.76 x 1O 9 M 
(observed mass M = 0.46 x 1O 9 M {5]). However, the BEC model gives a very good 
description of the inner region r < 1.5 kpc of the rotation curve. The NFW velocity fit fails 
to give any reasonable description of the rotation curve. 

• Holmberg II: Rather large velocity dispersions, as compared to what would be expected 
for a circular rotation velocity only, appear in the rotation curve of the dwarf galaxy Holm- 
berg II. Non-circular motions and strong baryonic gravitational perturbations in the galaxy 
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are the most probably explanation for the pattern of the rotation curve for r < 5 kpc. With 
the values of the radius and central density adopted to predict the rotation curve the mass 
of the galaxy is obtained as M§g?j = 2.76 x 1O 9 M (observed mass M = 2.07 x 1O 9 M 
jlj]). The BEC model describes well the outer (r > 5 kpc) region of the galaxy, but fails 
to describe the inner galactic regions, where any realistic model must take into account the 
presence of supplementary (local) interactions. The NFW fit cannot reproduce consistently 
the observed data. 

• DDO 53: The rotation curve of the galaxy DDO 53 shows a clear and smooth rotational 
pattern, even that the possibility of external perturbations for r > 2 kpc cannot be excluded. 
These perturbations may be responsible for the sudden decay of the curve, without the 
presence of a significant plateau. However, the BEC prediction describes well the rotation 
curve. The mass of the galaxy is obtained as M§qjj = 0.55 x 10 9 M Q , while the observed mass 
is M = 0.45 x 1O 9 M |5|. The NFW profile describes well the inner region of the rotation 
curves, but fails to give a good description of the curve in the outer galactic regions. 

• M81 dwB: The BEC model gives a very good description of the rotation curve of the 
galaxy. This galaxy is the smallest in the sample, having a mass of only M = 0.3 x 1O 9 M 
5j. With the adopted values for radius and central density the mass of the galaxy is 
M Mi?dwB = 0-33 x 1O 9 M . The velocity field shows no significant non-circular motions in 
the galaxy, except perhaps in the very outer regions. The NFW fit cannot reproduce the 
main features of this rotation curve. 

• DDO 39: The rotation curve of the galaxy DDO 39 shows the presence of significant 
non-circular motions in the inner galactic region with r < 1.8 kpc. A strong gravitational 
perturbation may be responsible for the distortion of the inner part of the rotation curve. 
Other non-circular motions may be present in the outer region of the curve, for r > 9 
kpc. The BEC model gives a relatively good description of the observational data in the 
middle and outer regions, where the motion is mostly rotational. The NFW fit gives a good 
description of the rotation curve. The adopted radius and central density for the BEC model 
gives a mass of M^g 39 = 1.6 x 1O 1O M . 

• DDO 154: The rotational velocity curve of the galaxy DDO 154 shows a regular rotation 
pattern. No significant non-circular motions can be observed in the galaxy. The rotation 
curve resembles that of a galaxy with solid-body rotation, with a steep increase in the inner 
regions. The BEC rotation curve is obtained by taking R = 7.5 kpc and p c = 6 x 10 -25 
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g/cm 3 , which gives a mass of M^q 154 = 4.71 x 10 9 M Q . The mass of the galaxy determined 
from observations is M = 5.40 x 1O 9 M js|. The BEC prediction describes reasonably well 
the outer region of the rotation curve. The NFW profile cannot give a satisfactory fit to the 
observational data. 



B. Density profiles and mean logarithmic density and velocity slopes 

The corresponding density profiles of the dark matter halos are presented in Fig. [2] for the 
BEC model, while the NFW density profiles are shown in Fig. [3J The density profile of the 
BEC condensate is non-singular, and well behaved at the center of the galaxy. Depending 
on the numerical values of the radius and central density, the mass distribution of the 
gravitationally bounded condensate decreases slowly as a function of r, with most of the 
matter concentrated in a core-like region. This theoretical prediction of the BEC model is 
in sharp contrast with the cusp-like structure found in the numerical simulations [2[, which 
can be seen in Fig. [3J The NFW density profile is singular at the center of the galaxy, 
for r — > 0, and it increases very rapidly for small r. On the other hand, for large r the 
NFW profile tends more rapidly to zero than the BEC profile. Another important difference 
between the two profiles is that the BEC profile predicts a finite and well defined radius for 
the dark matter distribution, whose density tends rigourously to zero on a surface that define 
the radius R of the dark matter halo. For the NFW density profile a well-defined radius does 
not exists. Therefore Bose-Einstein condensation processes in dwarf galaxies can affect the 
dark matter distribution in such a way that the central cusps predicted from dark- matter- 
only cosmological simulations are flattened, resulting in dark matter halos characterized by 
a large and almost constant density core, as observationally found in normal dwarf galaxies 
in the local universe 0, Q . 

As a next step in our analysis we consider the behavior of the logarithmic density profile 
slope a, defined as 

« = — , 15 

ctlogr 

and which can be easily expressed as 

7i 2 V 2 (r) 

a (r) = - [1 _ kr oot(fcr)] = (16) 
For r — >■ R, a(r) diverges, a (R) — > oo. The mean value of the logarithmic inner slopes of 
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the mass density profiles have been obtained from the observations of the rotations curves of 
dwarf galaxies p] and low surface brightness galaxies ((J, respectively. We define the mean 
value of the inner a as 



= -5- 



a(r)dr — — I (1 — ixx cot T(x) dx, 



(17) 
R in /R. 



where Ri n is the radius of the inner core of the dark matter distribution, and Xi 
The variation of (<Xj n ) as a function of is represented in Fig. HI 

For R in = 0.6R we obtain = —0.2818, a value which is very close to the value 

—0.29 ± 0.07 obtained in j^j. Therefore in the BEC dark matter model a value of (ai n ) in 
agreement with the observations can be easily obtained. Of course, the procedure proposed 
in the present paper for the determination of (aa n ) is very different from the method used in 
j^] , where first a break-radius, where the slope changes most rapidly, was determined. Then 
the inner density slope was measured by performing a least squares fit to the observational 
data points within the break-radius. In the theoretical estimation of we have assumed 
that Ri n = 0.6R is the break-radius. In the CDM model the value of a is a = — 1. 
The logarithmic slope of the rotation curves is defined as 

dhgV 



alogr 

and in the Bose-Einstein condensate dark halo model it can be obtained as 



(18) 



0(r) 



2^,2 



1 



k 2 r 



- 5 1- 



1 — kr cot(Ar) 
4:7rGpr 2 \ 



V 2 (r) 



(19) 



Therefore for Bose-Einstein condensed dark matter halos the logarithmic velocity and density 
slopes are related by the simple relation 



(5{r) 



1 + 



7T 



r y 



(20) 



a(r) \RJ 

This theoretical prediction represents another possibility for a direct observational test of 
the condensed dark matter model. 



IV. DISCUSSIONS AND FINAL REMARKS 



The Bose-Einstein condensation takes place when particles (bosons) become correlated 
with each other. This happens when their wavelengths overlap, that is, the thermal wave- 
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length Xt = y2iih 2 /rnksT is greater than the mean inter-particles distance a, At > a. 
The critical temperature for the condensation to take place is T cr < 2nh 2 n 2 ^ 3 /rnks 29 ]. 
On the other hand, cosmic evolution has the same temperature dependence, since in an 
adiabatic expansion process the density of a matter dominated Universe evolves as p oc T 3//2 
19| . Therefore, if the boson temperature is equal, for example, to the radiation tempera- 
ture at z = 1000, the critical temperature for the Bose-Einstein condensation is at present 
T cr = 0.0027-fT [19]. Since the matter temperature T m varies as T m oc a -2 , where a is the 
scale factor of the Universe, it follows that during an adiabatic evolution the ratio of the 
photon temperature T 7 and of the matter temperature evolves as T 7 /T m oc a. Using for 
the present day energy density of the Universe the value p cr = 9.44 x 10~ 30 g/cm 3 , BEC 
takes place provided that the boson mass satisfies the restriction m < 1.87 eV [l6|]. Thus, 
once the temperature T cr of the boson is less than the critical temperature, BEC can always 
take place at some moment during the cosmological evolution of the Universe. On the other 
hand, we expect that the Universe is always under critical temperature, if it is at the present 



time 



19]. 



Another cosmological bound on the mass of the condensate particle can be obtained as 
m < 2.696 (gd/ g) (Td/T cr ) 3 eV j^Q], where g is the number of internal degrees of freedom 
of the particle before decoupling, is the number of internal degrees of freedom of the 
particle at the decoupling, and is the decoupling temperature. In the Bose condensed 
case T d /T c < 1, and it follows that the BEC particle should be light, unless it decouples 
very early on, at high temperature and with a large g^. Therefore, depending on the relation 
between the critical and the decoupling temperatures, in order for a BEC light relic to act 



as cold dark matter, the decoupling scale must be higher than the electroweak scale 30 1. 

In the Bose-Einstein condensate dark matter model the physical properties of the galaxies 
are determined by three parameters: the scattering length a, the mass of the particle m, 
which completely determine the radius of the dark matter distribution R = R(a,m), and 
the central density of the dark matter, which is necessary for the estimation of the mass 
M = M (a, m, p c ). In both these quantities a and m enter in the combination a/m 3 . If the 
radius of the galactic dark matter distribution is known, then the ratio of the two intrinsic 
parameter that determine the properties of the Bose Einstein condensed dark matter halos 
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can be obtained as 

For the considered sample of dwarf galaxies the ratio a/m 3 varies between 3278.38 fm/eV 3 
(obtained from the galaxy M81 dwB with radius R = 1 kpc) and 504,084 fm/eV 3 for 
the galaxy IC 2574, with radius R = 12.4 kpc. Since the variation of the radius for the 
considered sample of galaxies is around one order of magnitude, the range of a/m 3 obtained 
from individual galaxies can differ by two orders of magnitude. By taking the mean value of 
the radius Rmean = 6.5 kpc, we obtain a/m 3 = 122189 fm/eV 3 . A very small particle mass, 
of the order of 10 -22 eV would require a scattering length of the order of a ~ 1.22 x 10~ 61 
fm, which seems to be too small to be realistic. The variation of the mass of the dark matter 
particle as a function of the scattering length a is represented, for each considered galaxy, 
in Fig. [5j 

The contour plot of the mass of the dark matter particle as a function of the galactic 
radius R and of the scattering length a, m = m (a, R) is represented in Fig. |6j A contour 
plot, which is an alternative to a three dimensional surface plot, is a graphical technique 
for representing a three-dimensional surface by plotting constant m slices, called contours, 
on a two-dimensional format. That is, for a given value of m, the curves shown in the 
figure connect the points with coordinates (R, a) where that m value occurs. In order 
to consistently derive the values of a and m from observations another relation between 
these two parameters is necessary. Such a relation could be obtained from cosmological 
considerations, as suggested in [281 ] . 

Finally, we would like to address the problem of the mass of the dwarf galaxies. As 
compared to the standard estimation of the mass M w AnR 3 p c /3, the total mass of the 
Bose-Einstein condensate, given by Eq. (TTTj) . is smaller by a factor of three for constant 
density halos, and can be one order of magnitude smaller for halos with a non-constant 
density profile. This may solve the mass discrepancy between observations and simulations, 
which was pointed out in [7J]. As compared to the BEC models, the usually core-like models, 
such as pseudo-isothermal halo models, dominated by a central const ant- density core jlj], 
predict a much higher galaxy mass. 

Since dwarf galaxies are usually considered to be dark matter dominated, the possibility 
of reproducing their rotation curves, and the mean value of the logarithmic inner slope of 
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the mass density profiles, represents a significant evidence in favor of the strongly - coupled 
dilute Bose-Einstein condensate dark matter model. Moreover, the density profile of the 
condensate has generally a large core with a slowly decreasing density. Another important 
advantage of the model is that all physical quantities are well behaved for all r, and there 
are no singularities in any of the physical parameters. In the present paper we have shown 
that the strongly correlated Bose-Einstein condensate dark matter model gives generally a 
relatively good description of the properties of the dark matter dominated dwarf galaxies. 
All the relevant physical quantities can be predicted from the model, and can be directly 
compared with the corresponding observational parameters (the dark halo mass, the radius 
of the galaxy, and the observed flat rotational velocity curves, mean values of the density 
profiles etc.). Therefore this opens the possibility of in depth testing of the Bose-Einstein 
condensation models by using astronomical and astrophysical observations at the galactic- 
intergalactic scale. In this paper we have provided some basic theoretical tools necessary for 
the comparison of the predictions of the condensate model and of the astronomical results. 
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FIG. 1: Comparison of the velocity of the rotat[i§n curves predicted by the Bose-Einstein Conden- 
sate dark matter model (solid curve), of the velocity fits obtained with the NFW density profile 
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FIG. 2: Density profiles of the galaxies IC 2574 (solid curve), NGC 2366 (dotted curve), Ho II (short 
dashed curve), DDO 53 (dashed curve), DDO 39 (long dashed curve) and DDO 154 (ultra-long 
dashed curve), respectively. 
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FIG. 3: NFW density profiles of the galaxies IC 2574 (solid curve), NGC 2366 (dotted curve), 
Ho II (short dashed curve), DDO 53 (dashed curve), DDO 39 (long dashed curve) and DDO 154 
(ultra- long dashed curve), respectively. 
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FIG. 4: The mean value of the logarithmic inner slope of the mass density profile of dwarf galaxies. 
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FIG. 5: Mass (in a logarithmic scale) of the dark matter particle as a function of the scattering 
length a for the considered sample of galaxies, in increasing radius order: M81 dwB (solid curve), 
Ho I, DDO 53, NGC 2366, DDO 154, DDO 39, and IC 2574, respectively. 
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FIG. 6: Contour plot of the mass m of the dark matter particle as a function of the scattering 
length a and of the radius R of the dwarf galaxy. 
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